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A set of quasi-exactly solvable quantum mechanical potentials associated with the Poschl- Teller 
potential, the generalized Poschl- Teller potential, the Scarf potential, and the harmonic oscillator 
potential have been studied. Solutions of the Schrodinger equation for each potential have been 
determined and the eigenstates are expressed in terms of the orthogonal polynomials. The potentials 
' are related to each other by suitable change of variables. 
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INTRODUCTION 



The discovery of new class of physically significant spectralproblems, called quasi-exactly solvable(QES) models^, 
2;], has attracted much attention^ U Several methods^ |(| for the generation of the QES potentials have been 
worked out and consequently a number of QES potentials have been suggested. One of the methods for the calculation 
of eigenstates and eigenvalues of the QES potentials is the use of orthogonal polynomials. Bender and Dunne have 
showed0 that there is a correspondence between the QES models in quantum mechanics and the set of the orthogonal 
polynomials P m (E), which are polynomials in energy E. In particular they have demonstrated the properties of the 
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polynomials for the QES sextic oscillator potential[8|. Other aspects of the P m (E) have been discussed in various 
' papers 0, In this paper we show that the polynomial P m {E) normally appears as a factor of the coefficient 

in the expansion of the eigenstates of the Schrodinger equation in the presence of the QES Poschl- Teller potential, 
the Generalized QES Poschl- Teller potential and the PT symmetric Scarf potential. We also show that the P m (E) 
O , satisfies the same three-term recurrence relations for all three potentials. 

An algorithm generating the analytical eigenfunctions as well as the eigenvalues of the Schrodinger equation for 
various QES potentials is developed. The procedure presented here reproduces the results of the exactly solvable 
Schrodinger equations in proper limits. Our formalism s igna ls that, there is a closer similarity between QES Poschl- 
Teller potential and the perturbed Poschl- Teller potential |lllll2| . It is well known that the harmonic oscillator and the 
Poschl- Teller potential possess the similar theoretical behavior in a certain range of variable. In contrast to enormous 
interest in the perturbed harmonic oscillator, there have been only a number of studies on the perturbed Poschl- Teller 
potential [Tl| 

We show that the QES Poschl- Teller potential can be transformed on to the generalized QES Poschl- Teller potential 
and the QES PT symmetric Scarf potential by replacing coordinate x — > x/2 and x — > a;/2 + «7r/4, respectively. These 
potentials are also related to the QES harmonic oscillator potential by redefining the parameters and applying an 
appropriate limiting procedure. 

The paper is organized as follows. In section 2 we discuss the generation and solution of the QES Poschl- Teller 
potential. The QES generalized Poschl- Teller potential and the QES PT symmetric Scarf potentials are obtained from 
the QES Poschl- Teller potential in sections 3 and 4, respectively. Transformations of QES Poschl- Teller potential to 
the sextic oscillator potential and the QES PT symmetric Scarf potential to the PT symmetric anharmonic oscillator 
potential are discussed in sections 5 and 6 respectively. In conclusive remarks we discuss the implementation of our 
method on the other QES potentials which will be the topic of another publication. 



QES POSCHL-TELLER POTENTIAL 



The QES Poschl- Teller potential can be generated by several methodspil fill Il3|. One method is to use the Lie 
algebraic technique. The linear and bilinear combinations of the operators of the sl(2,R)^^| Lie algebra with the 
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standard realization|i|,0,^|], leads to the following differential equation, 

*(1 - ^^r 1 + (L + l + z(B + 4j - lA 2 z)) d ^- - (A - 2jgA 2 z)5R,(z) = (1) 

where L, g, A and A are constants and j = 0, 1/2, 1, • • • . The differential equation Q becomes QES, provided B is 
taken as: 

B = -- (2L + 8j + 5 - y/l + AA{A + \ + {2L + 8j + 5)gA)) . (2) 

Moreover, it becomes exactly solvable with the condition q = 0. The function ^Stj(z) is a polynomial of degree 2j. In 
order to obtain the QES quantum mechanical potentials we can transform Q in the form of Schrodinger equation by 
introducing the variable 

z = — sinh 2 ax. (3) 

Now we define the wave function as follows: 

i/)(x) = (coshax) qA2 - B - L - i: >- 2 (smhax) 1+L e-^ A2 cosh2Q2 ^(- sinh 2 ax). (4) 
Substituting © and Q in we obtain the Schrodinger equation with (H = 2m = 1): 

d 2 i){x) 



dx 2 

where the potential V(x) is given by 



+ {V(x) - E)^(x) = (5) 



V(x) = L(L + l)a 2 esc h 2 ax - A(A + l)a 2 sec h 2 ax + 

q(2BA 2 sinh 2 ax + qA 4 sinh 4 ax)a 2 tanh 2 ax (6) 

The eigenvalues of energy are given by the expression 

E = [-(L -A + 2m + l) 2 + (L - A + B + Aj + 2)(2L + 4m + 3) + 4A] a 2 . (7) 

One can check that for q — the potential given in © is exactly solvable Poschl- Teller potential and the eigenstates 
of the Schrodinger equation can be expressed in terms of the Jacobi polynomials. It is easy to see that the equalities 
B = A — L — 4j — 2 and A = im(L — M + m + 1) hold when q = 0. Then the eigenvalues of the Schrodinger equation 
take the form 

E = -a 2 [{L - A + 2m + l) 2 + 4m(L -M + m+ 1)] . (8) 

Next task is now to determine the eigenfunction of the QES Schrodinger equation(|!)J|. Therefore, we search for a 
solution of by substituting the polynomial 

2j 

5%(z) = ]T a m z m (9) 

m=0 

which leads to an expression for the coefficients a m 

_ (4 g A 2 r(2j)!(2L+l)!(L + m)! D 
m ~ 2m!(2j-m)!(2£ + l + 2m)! m( h 1 1 

The polynomial P m (A) satisfies the following three-term recurrence relation 

2(2j - m)gA 2 P m+ i(A) + m(2L + 2m + l)P m _i(A) - 2(A + m{B + Aj-m+ l))P m {\) = (11) 

with the initial condition Pq(X) — 1. The polynomial P m (A) vanishes for m ^ 2 j + 1 and the roots of P2j+i(A) = 
corresponds to the A— eigenvalues of the Schrodinger equation JSJ. The first three of these polynomials are given by 

Pi = A 

P 2 = X 2 -(B + 4j)X-j{2L + 3)qA 2 

P 3 = A 3 — (3P12j - 4)A 2 + (12) 
[2B(B - 2) + 16j(B + 2j - 1) + {2L - j(6L + 13) + 5)qA 2 ] X + 
2j(2L + 3){B + 4j -2)qA 2 
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The recurrence relation (|ll|l can also be put in the matrix form. The tridiagonal matrix characterizes the system, 
//3o-4A n aj \ ( P \ 



7i 



Pi - 4A fi 2 



Po 
Pi 



= 



(13) 



\ 72, 0y-4A/ V 7 

where the parameters in matrix elements are given by 

j m = 2m(2L + 2m+ 1), /i m = 4mqA 2 , and (3 m = 4m(B + 4j - m + 1). (14) 

Analytical solutions of the recurrence relation (|11|) and the determinant of l|13|) for A are available only for the first 
few values of j ^ 2. For j > 2 the solutions become numerical, and the numerical errors grow rapidly. The solutions 
take simpler forms and the precision becomes better for A ^> q, in which case B takes the value 

B s» 2(A + m(L - A + m + 1)) - m(2L + 8j + 5)gA (15) 

In this approximation the QES Poschl- Teller potential in Q becomes comparable with the perturbed Poschl- Teller 
potential, a specific form of which has been studied in [15| . 



GENERALIZED QES POSCHL-TELLER POTENTIAL 

In this section we present a procedure that relates the QES Poschl- Teller potential to the generalized QES Poschl- 
Teller potential. It is amusing to observe that when the coordinate x of the Poschl- Teller potential is replaced by 
x — > x/2, the QES Poschl- Teller potential transforms to the generalized QES Poschl- Teller potential: 

2 

V(x) = y[(L(i+ 1) + A(A + l))csch 2 ax 

+ (L - A) (L + A+ l)coth ax csc hax] (16) 
+q{2BA 2 smh 2 — + qA 4 smh 4 -y)" tanh^ — 

Obviously when q — this potential reduces to the exactly solvable generalized Poschl- Teller potential. The wave 
function corresponding to the solution of (J5J with the potential of l|16|) now takes the form 

^(s) = (cosh^)« A2 - B - L - 4 -''- 2 (sinh^) 1 + L e-^ A2cosh « 2: 3?,(-sinh 2 ^). (17) 

The corresponding energies are determined as 

E 

E' = — + m(B + ij — m+ l)a 2 (18) 

where E is the eigenvalue of the Schrodinger equation with the QES Poschl- Teller potential given byQ. It is clear 
that the same recurrence relation in (|ll|) for the polynomial P m (A) holds true. 



QES PT SYMMETRIC SCARF POTENTIAL 

It is interesting to observe that one can transform the QES Poschle- Teller potential to the PT symmetric QES Scarf 
potential by replacing x by the complex variable 

x in , 
X ^2+- a ^ 

Then one obtains the following potential: 

2 

V(x) = - y [i L i L + 1) + MA + 1)) sec h 2 ax 

+ i(L - A) (L + A + l)tanh ax sec hax] (20) 

a 2 .„ _ .9 . , o,ax Z7T. 9 ,a . ,A,ax i7r.. , <>,a.x in, 
+ -jVdBA 2 sinh 2 (— + _) + (? 2 A 4 smh 4 (— + _))tanh 2 (— + — ) 



4 



This is the QES form of the PT symmetric Scarf potential which has not been discussed in the literature. The exactly 
solvable part of the potential H20[) fgiven in the [...] which, can be obtained if q = 0) has been discussed recently by 
Bagchi 14]. The energy eigenvalues are the same as in (|18|) obtained for the generalized QES Poschle- Teller potential. 
But the wave function now reads 



/»+tanhrvT\ 5( B +2£-<?^ 2 +4j+3) 



e-i« A2sinhax 5i i (-sinh(— + ^). (21) 

Here again the polynomial P m (A) satisfies the same recurrence relation defined in llll) 

The method described in this section in order to transform one type of potential to another type may be further 
generalized. The transformation x — > ax+ib , where a and 6, are arbitrary real parameters, preserves the PT symmetry. 
The energy eigenvalues of the new potential involves further terms in addition to the scaled energy eigenvalues of 
the QES Poschl- Teller potential. The corresponding wave function is obtained from the former wave function of the 
Poschl- Teller potential by replacing x — > ax + ib, and the polynomial P m (A) obeys the same recurrence relation of (|llll . 
This transformation has been discussed in the literature 15] for the exactly solvable potentials. 

THE SEXTIC OSCILLATOR 

In this section we discuss a method about the transformation of the QES Poschle- Teller potential to the sextic 
oscillator potential. The QES Poschle- Teller potential can be converted to the radial sextic oscillator potential by 
redefining the parameters and taking suitable limits while keeping the variable x intact. If we redefine the parameters 
in JBJ and @ by introducing 

b , a e — mb , . 

B=—, A=—, A = — (22) 

a a a 

and taking the limit of the {V(x) — E) term in the Schrodinger equation for a — ► the Poschle- Teller potential 
transforms to the radial sextic harmonic oscillator potential: 

V = Li - L t 1) + (b 2 - {2L + 8j + h)qa 2 )x 2 + 2bqa 2 x 4 + q 2 a 4 x 6 . (23) 
x z 

The ground state wave function^) takes the form 

iP(x) = x 1+L e-i x2 -i x * (24) 

and the energy eigenvalue is given by 

E = Ae + (2L + 3)6. (25) 
The general solutions of the Schrodinger equation with the sextic oscillator potential can be obtained by substituting 

ip(x) ^x 1+L e-i x2 -^^ J {x 2 ) (26) 
into the Schrodinger equationJSJl. By rearranging the terms one can obtain, 

n ; ^ 2m\(2j - m)\(2L + 1 + 2my ' 



m=0 



along with the recurrence relation satisfied by P m (e) 

2(2j - m)qa 2 P. m+1 (e) + 2(e - bm)P m (e) - m{2L + 2m + l)P m _i(e) = 0. (28) 
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The first four polynomials are given by 
Pi(e)=e 

P 2 (s) = e 2 - be - ~(3 + 2L)qa 2 

P 3 (e) = e 3 - 3be 2 - 2{b 2 - 2{L + 2))e + 26(3 + 2L)qa 2 (29) 
Pi{e) =e 4 - 6be 3 + (lib 2 - 5(2L + 5)qa 2 e 2 + 3b(-2b 2 + (10L + 2l)qa 2 e 

~b 2 (2L + 3) + i(4L(5L + 1) + 2\)q 2 a A 

They agree with the polynomials given in references and [13. We have also compared our results with the 
numerical solutions obtained in the reference 16], for the potentials, 

Vi{x) = x 2 + 



2(7.625) 3 / 2 7442 

^^ 2+ 2(7^ + ^2 <™ 

x 4 

Vsix) = ^ + 2(7.125) 3 / 2 + 6498 
where we have obtained exactly the same results of |l6j . for the energy eigenvalues 

E x = 2.897143, E 2 = 5.891677, and E 3 = 8.991223. (31) 
corresponding to the potentials V\{x), V 2 {x) 1 and T/3(x)respectively. 

PT SYMMETRIC ANHARMONIC OSCILLATOR POTENTIAL 

In order to transform the potential given in (|2l)ll to the PT symmetric anharmonic oscillator potential we replace 
q and L by 

9->^> L^±(qL 2 + (2 + 4A)L-(20-32j)). (32) 

Further, substituting in the parameters in l|18fl and (|20|l 

, — qa 2 7 b If , — / 7 \ ab 



.4 



a:' 



3 ■frfa 2 ' ^ 2 \ \ y/Vf J qa 2 



e + 2jb ijqa 2 4qa 2 If 7 \ 2b 2 + 17iqa 2 
a z or or 17 \ v 17 / qa z a 

and taking the limit of the (V(x) — E) given in l|18|) and (|20|l when a — > 0, we obtain the PT symmetric anharmonic 
oscillator potential: 

V = 2i(b£ - (1 + 2j)qa 2 )x + {b 2 - 2lqa 2 )x 2 + 2iqba 2 x 3 - q 2 a 4 x 4 . (34) 
This leads to the energy eigenvalues, 

E = e + b(l + 2j)+P. (35) 

The ground state wave function of the potential can be obtained from l|21|) by using the same limiting procedure as 
introduced in (|32"|) and lj3lfl) : 

$ = e (—^ 2 -¥ ^ 2 * 3 ) ( 36 ) 
The wave function for any j can be obtained by letting 

i,( x ) = e (- lx -kb* 2 -% q A^ 3 ) a m x 2m . (37) 
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Here we obtain a four-term recurrence relation for the energy, 

2i(2j - m)qa 2 P m+1 (e) + (e - 2b(m - j))P m (e) - 2imiP m -i[e) + m(m - l)P TO -a(e) = (38) 
The first four P m is given by 

Pi(e) = e 

P 2 (e) = e 2 -b 2 -Aqa 2 £ 

P 3 (e) = e 3 - 4(b 2 + \qa 2 ()e - 16gV (39) 
P 4 (e) = e i -10{b 2 + 4qa 2 l)e 2 -96q 2 a 4 e + 9{b i + 8qb 2 a 2 i+16q 2 a 4 e 2 ) 

The results are in agreement with those given in|l7|. 



CONCLUSIONS 



We have developed a general procedure to obtain the eigenstates of some QES potentials in terms of the orthogonal 
polynomials as well as the eigenvalues. We have proven that these potentials can be obtained from each other by 
adequate transformations. Some examples have been presented to test the validity of the procedure given here. 

The present work can be generalized in various directions. The method we have introduced can be applied to 
determine the QES forms of the other exactly solvable potentials. The present construction seems to exhibit closer 
similarity to the current perturbation constructions. In particular, the perturbed Poschl- Teller potential can be 
worked out by following the procedure discussed here and some of its eigenstates can be determined exactly. 
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